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Abstract 

We study the diffusion of a linear polymer in the presence of permeable membranes 
without excluded volume interactions, using scaling theory and Monte Carlo simulations. 
We find that the average time it takes for a chain with polymerization index N to cross 
a single isolated membrane varies with N as N^^"^, giving its permeability proportional 
to N'^ . When the membranes are stacked with uniform spacing d in the unit of the 
monomer size, the dynamics of a polymer is shown to have three different regimes. In the 
limit of small d <^ N^^^, the chain diffuses through reptation and D ^ N^^. When d is 
comparable to N-^^^ the diffusion coefficients parallel and perpendicular to the membranes 
become different from each other. While the diffusion becomes Rouse-like, i.e. D ^ N^^, 
in the parallel direction, the motion in the perpendicular direction is still hindered by the 
two-dimensional networks. The diffusion eventually becomes isotropic and Rouse-like for 
large > iV. 



1 Introduction 

The physics of tethered membranes |2|, |^, ^ has been studied with great interest due not only 
to its relevance to biology but also to its numerous applications in technology. Stabilizing 
liposomes, which are used in drug delivery, by tethered membranes is just one of many 
important applications B. It is also possible that tethered membranes could provide a new 
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way of filtrating polymers in a solution according to their molecular weights. This method 
should be especially useful where conventional methods such as gel electrophoresis ^, |8| do 
not work. 

Motivated by this, we study in this paper the effects of permeable membranes on the dy- 
namics of a flexible linear polymer. More specifically we focus on the diffusion characteristics 
of a linear polymer in the presence of a single or multiple membranes. 

We model our permeable membranes as two-dimensional regular networks of polymers 
with the mesh size bigger than the persistence length of the diffusing chain. This type 
of membranes can be made by crosslinking linear polymers using ultraviolet light |^ or by 
extracting spectrin networks from human red blood cells Q. The widely used molecular filters, 
called zeolites [10], may also provide a good system if the molecular parameters such as the 
Kuhn length of the diffusing polymer is sufficiently small. 

For simplicity, we restrict our discussion to systems with no excluded volume interactions. 
The only constraints the two-dimensional networks impose on the dynamics of the chain 
polymer are topological barriers. From a theoretical point of view, this is a new problem quite 
distinct from the study of reptating chains in an inhomogeneous environment [pT|, which 
has been investigated extensively partly due to its relevance to DNA gel electrophoresis. Even 
though the main ingredients in reptation theory are the topological constraints imposed by 
surrounding chains or gel networks, most studies on reptation in an inhomogeneous medium 
concentrates on more complex (and presumably more realistic) systems with excluded volume 
interactions. 

The remainder of this paper is organized as follows. In section ^ we first define our 
problem in more precise terms and give a brief account of the simulation methods used in 
this work. In section ^, we consider the longest relaxation time of a linear polymer in the 
presence of a single membrane. This is obtained numerically by measuring the first passage 
times of the linear polymer through an array of membranes which are placed with a spacing 
proportional to the radius of gyration of the polymer. In section ^, we use the data obtained 
in the previous section to calculate the permeability for our two-dimensional network, which 
is one of the most important characteristics of permeable membranes. In section |^, we study 
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the multiple membrane case with uniform spacings in more detail. The diffusion constants 
parallel to the membranes as well as the perpendicular ones are studied. We first present a 
simple scaling theory as a function of the polymerization index and the spacing between the 
membranes. Then we present the simulation results which support this analysis. In section^, 
we consider the mean square displacements of a monomer as a function of time when a chain 
disengages from a network. Finally in section |^ we summarize our results and discuss possible 
applications of this study. 



2 The Model 

We consider a random walk chain of N monomers or beads on a cubic lattice with a lattice 
spacing 6, which is the monomer size or the bond length of the chain. We only consider 
even A^'s due to the bipartite nature of the cubic lattice. 

The simulation is done by allowing the chain to move according to the following set of 
Monte Carlo rules which mimics Rouse dynamics. 

I Two consecutive bonds which form an L-shape can be inverted, in the plane formed by 
the L. 

II A kink, two consecutive bonds which are on top of each other, can be flipped to one 
of the six lattice directions including the original one. (See, for example, Fig. 1 of 



reference 1131.) The free ends of the chain are treated as kinks. 



Now we consider another cubic lattice with lattice spacing a which interpenetrates the 
first one. We assume a is an integral multiple of b so that the lattice points of the second 
lattice always reside at the centers of the unit cells of the first. We let a = 6 = 1 in the 
numerical simulations. 

Our membrane is a two-dimensional slice of the second lattice in the horizontal plane. 
Since the membrane occupies a different lattice than that of the chain, it does not alter the 
statics of the chain. But we assume that the chain and the membranes cannot pass through 
each other. That is, the links of the membranes act as topological barriers for the motion of 
the chain. 
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This is implemented in our simulation by modifying rule I. We allow the L-inversion only 
when the chain does not pass through a link of a membrane by this move. This is basically 



the same as the Doi-Evans-Edwards algorithm|14, 15|. Note that the dynamics of a chain is 
altered only near the vicinity of the membranes and if the plane formed by the L is horizontal 
rule I is not affected at all because the links of the membranes are all horizontally placed. 



3 Relaxation time of chain penetration 

First we study the effects of a single membrane on the dynamics of a linear polymer with 
polymerization index A^. One of the most basic questions we have to ask in this situation is 
how much the membrane will slow down the diffusion of the polymer which is in its vicinity. 
More specifically we ask what is the time scale of the complete penetration of the polymer 
through the two-dimensional network. 

Before we define the relaxation time for this penetration, we point out some problems 
associated with the very concept of the "complete penetration" . First we note that because 
a polymer has a finite size of order the radius of gyration, Rg ^ N^^'^b, the presence of the 
membrane interferes with the Rouse dynamics of the chain as long as its center of mass 
lies in a slab-like region around the membrane. We call this region of thickness ~ -Rg the 
entanglement zone. Here we stress that the thickness of the entanglement zone is a function 
of N. 

If we try to define the relaxation time of a chain through a membrane by diffusion, we 
soon encounter difficulties. Diffusion is a random process and the polymer may never pass 
through the membrane. If the polymer penetrates the membrane at some point in time, the 
question remains as to how much time should be attributed to the effect of the membrane. 

We circumvent this problem in the following way. First we consider two configurations, 
one with the center of mass of the polymer at a distance ~ Rg from the membrane and 
another in which the polymer is maximally entangled (on average) with the membrane, that 
is, its center of mass lies at the membrane. Since the equilibrium statics are not affected by 
the topological membrane, these two configurations should arise with the same probability 
in equilibrium. Then from the principle of detailed balance, the transition probability from 
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one configuration to the other should be the same as the reverse. Hence the transition times 
in both directions are the same. 

If we place a chain in the maximally entangled position (as defined in the previous para- 
graph) the chain eventually diffuses out to either end of the entanglement zone. We define 
the average first passage time in this process as the escape time. Then we define the pen- 
etration time as (4 times) the escape time since, as described above, the escaping process 
takes the same amount time as the reverse process and the penetration is, by definition, the 
combination of the two. Note that the chain does not necessarily penetrate to the other side 
of the membrane in this time scale. 

The penetration time defined this way can be estimated easily by simple application 
of the reptation theory jl^. First we note that the chain typically crosses the membrane 

iV^/^ times in the entanglement zone with the average contour length ~ N^^'^b between 
the crossing points, or the entanglement points. The relaxation of the chain can be divided 
into two processes; Rouse relaxation of chain segments between entanglement points, which 
we call Rouse segments, and diffusion of kinks past entanglement points. 

Even though this situation is quite different from the one of the usual reptation theory 
described by "reptilian" motion in a tube, we can still imagine the chain wiggling through 
a series of entanglement points inside an imaginary tube, its mean diameter being of or- 
der N^/^b. 

The longest relaxation time in reptation theory is tiN^/Nq, where A^e is the number of 
monomers between entanglement points. Since ~ N^^'^ in our case, we obtain|lf:] 



iv^/Vi (1) 



Therefore dynamics of a chain through a tethered membrane is faster than the reptation in 
a three-dimensional network Trep ~ N^ti and slower than the Rouse diffusion tr, ~ N'^ti, as 
expected. Note that A^^/^-dependence has been known to exit in various contexts such as the 



dynamics of a ring polymer in a gel 1 17]. 



In a time scale of the chain moves a distance of the order Rg N^/'^b. Hence if we 
place an infinite array of membranes with ~ Rg apart from one another then the average 
time for a chain to hop across a membrane should be proportional to the relaxation time t^. 



We conducted a simulation under these conditions for chain lengths 8 through 80. The inter- 
membrane distances are set to 2bN^^'^, with 6 = 1. At these distances the chains touching 
two adjacent membranes at the same time is highly unlikely. The penetration time, t ~ Tm 
is numerically defined as follows. 

where Tp is the first passage time from the center of one inter-membrane zone to a neighboring 
one. The angular brackets indicate the average over the distribution of Tp's. 

The log-log plot of r vs. is shown in Fig. ||. The typical number of Tp's used in obtaining 
each point is approximately 10'^. As the figure shows the exponent is clearly close to 2.5 as 
the above scaling theory predicts. If there were no fluctuations in the size of the chain, we 
should be able to produce more accurate results by shrinking the intermembrane distance 
down to the radius of gyration of the chain. Because of the more generous spacing used in 
the simulation the results must contain a part from Rouse dynamics. But the effect should 
be small because » tr for large N. 

The penetration time is also expected to depend on the mesh size of the membrane, a. 
Even though the a dependence is harder to ascertain in our lattice models since a only takes 
integer values, it can be obtained using simple scaling argument as follows. First we assume 
Tjn is a function of only N and a/b. And since we already know the A-dependence we write 
Tin as 

Tn. = nA^/2 (^)' (3) 

where p is the scaling exponent. When a approaches N^^'^b the dynamics should be Rouse- 
like. Hence we get p = —1. That is, is linearly proportional to the inverse of the mesh 
size a unlike r^ep, which has the quadratic dependence on a~^. 



4 Permeability of a membrane 

As suggested earlier, two-dimensional networks of polymers or two-dimensional slabs with 
pores could be used for filtration purposes. One of the most fundamental quantities in the 
membrane filtration technique is the permeability of a membrane, P, which is defined by the 
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following empirical equation, 

J = -PAc (4) 

where J is the flux density of the solute molecules through the membrane and Ac is their 
concentration difference across the membrane. 

A membrane is called selectively permeable when P depends on the solute being filtered. 
Our model membrane is selectively permeable to polymer solutes of different lengths as will 
be shown below and this is the basis for our proposal that tethered membranes be used as a 
dialysis tool. 

From Fick's law we can derive P in terms of quantities which are more easily measurable, 
the diffusion constants. 

' ' (5) 



P J entanglement zone 

where is a diffusion constant in the perpendicular direction to the membrane and the 

thickness of the entanglement zone is of order Rg, as described earlier. The diffusion constant 

in this entanglement zone should vary due to the spherical shape of the polymer, as indicated 

by the z-dependent in Eq. (^). But we neglect this as in the previous section and define 

the average diffusion constant along the z-direction, D^, using the following formula 

-R(T f dz , , 

(6) 



J D,{z) 

Then from the definition of Tm we get 
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Hence from Eq. (^) or 



we get 



= Ri/T^ ~ N—2 (7) 



(9) 



5 Diffusion constants of a polymer 

Now we consider the multiple membrane case. We stack membranes along the z-direction 
with a spacing d. A linear polymer of polymerization index N is introduced at an arbitrary 
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position in this system and the diffusion constants are measured both along the z-direction 
and in the xy-plane using the fohowing equations. 

and 

where X, Y and Z are center of mass coordinates of the chain. 

First we expect that there should be a crossover in the D's at dc ^ Rg- When d is smaller 
than dc the chain always entangles with multiple membranes and we expect that some sort 
of modified reptation should occur. On the other hand when d is much larger than dc the 
membranes act as small perturbations to the chain dynamics which should be otherwise 
Rouse-like. We study both limits in turn. 

5.1 (i<4 

In this dense regime of membranes we expect the reptation theory to be useful in describing 
the dynamics of a linear polymer, since the array of closely placed membranes should not be 
much different from a three-dimensional network. 

We first estimate the density of entanglement points. The polymer intersects of order 
Rg/d ~ N^/'^b/d membranes on average. Furthermore it crosses each membrane of order 
Rg/a ~ N'^/'^h/a times, which can be seen by rescaling the monomer size h to the mesh 
size a, so that Nh"^ /a?. Hence the total number of entanglement points is R^/ad ^ 

Nb^/ad. Then the average number of monomers between entanglement points is obtained 
by dividing N by this number, that is, we obtain N^. ~ ad/lP'. Naive application of reptation 
theory 1 16 1 as in section H], yields 
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T = n-^, (12) 

ad 

D = (13) 



The diffusion constant D increases as a linear function of the spacing d. 
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The simulation results for D^, and D^y for = 50 are plotted in Fig. |2| against d up 
to d = 30. Each point in this figure was averaged over 1000 Monte Carlo runs, each spanning 
longer than the penetration time. As the above argument based on the reptation theory 
predicts, diffusion should be isotropic in this regime. As Fig. |2| shows, this is indeed the case 
in the limit of d/Rg 0. But the system becomes anisotropic rather rapidly as the spacing 
increases. While remains linear in d even well above dc, D^y deviates dramatically from 
a straight line around a relatively small d. 

The scaling behavior in this regime is shown in Fig. ^ and Fig. ^ for chain lengths up 
to 80. As d increase beyond dc the scaling clearly breaks down, which is the topic of the next 
subsection. 

5.2 (i ~ 4 

In this regime the chain diffuses in a layered structure of free space and membrane. Thus 
depending on the relative importance of each, the dynamics of the chain can be either more 
Rouse-like or more reptation-like. 

We can use a simple resistor-network analogy] 18] to derive the effective diffusion constants 



in this regime. As will be shown shortly and D^y have very different behavior. 

First we start with D^- We divide the space into two categories; an entanglement zone of 
thickness Rg and a Rouse zone of thickness d — Rg. The effective resistance per length d in 
the z-direction per unit (xy-) area, TZ, is the sum of two components, 

71 = 7^E + 7^R (14) 

or 

d R^ d — Rff ,^ 

where is the Rouse diffusion constant and ~ N^^/"^ is the average diffusion constant 
in the entanglement zone as obtained in earlier sections. Then we obtain 
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Note that the crossover point from the entangled to the Rouse behavior is 

d^R^^r^Nb (17) 

contrary to the naive expectation of N^^'^b. 

The diffusion constant in the xy-plane, D^y, can be obtained in a similar manner except 
that the two components divided as above are connected in parallel rather than in series. 

The effective conductance per unit length along a direction perpendicular to the z-axis, 
per rectangular area with height d and unit width, denoted by C, is again the sum of two 
components, 

C = Ce + Cr (18) 

or 

dD^y = RgDE + id- Rg)DR (19) 



yielding 



D,,, = D„|l-^(^i!__M) (20) 



The entangled dynamics cross over to Rouse dynamics at d ~ i?g differing from the case 
of above. Hence in the spacing range N^^H <d<Nb, we have very anisotropic diffusion 
constants. While the dynamics in the j;y-plane is almost free the effects of the membranes 
are still strong along the z-direction. Beyond this spacing ~ Nb, the dynamics becomes 
isotropic again with Rouse-like diffusion constants in both directions. The trend is clear from 
Fig. |. 

The permeability for a stack of m membranes can be defined in a similar way as in 
section The permeability is the diffusion constant perpendicular to the slab divided by 
its effective thickness as shown in Eq. (|8|). Since the thickness of the stack of m membranes 
is md for large m, the permeability P is given by 

For large m, P varies as N~'^ in the dense regime and in the dilute regime identical to 
the scaling of Dz- 
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6 Escape of an entangled chain from a membrane 



So far we have concentrated on the diffusion of the center of mass of a chain. Now we look 
at the mean square displacements of a monomer in the middle of the chain, which reveals 



small-scale dynamics in more detail [ 19 1. 

First we place a chain and a membrane so that the chain is maximally entangled, that 
is, the center of mass of the chain lies in the membrane. After the escape time (defined in 
section ^), the chain will diffuse from the entanglement zone into free space. This process is 
shown in Fig. ^ for N = 200. First we see t^^^ behavior up to Nti which should be the local 
relaxation of Rouse segments, that is, the chain segment relaxation between entanglement 
points. At greater times the membrane starts to restrict the motions of chain which makes 
the relaxation much slower, probably proportional to t^^^ as shown in the figure. In this 
regime the dynamics in the z-direction and xy-direction are different as expected because of 
the two-dimensional nature of the membrane. We need to simulate larger chains to obtain 
the fine structure more accurately in this regime. 

After the Rouse time scale N'^ti, the chain starts to escape from the membrane. This 
process is unusually fast in the sense that the exponent is larger than 1, about 3/2 as shown 
in Fig. ^. This is probably because the farther the chain moves from the membrane the less 
entangled it is with the membrane, and hence the larger the diffusion constant. The chain 
will eventually escape from the membrane at a time scale N^^'^ti as obtained in section ^ 
and normal diffusion (Z^) ~ t will resume, which is not shown in the figure. 

The monomer relaxation behavior for densely packed multiple membrane case can be 
interpolated from the data shown here and those given by reptation theory in a three- 



dimensional network. See, for example, reference |19 



7 Conclusion 

The dynamics of a flexible polymer in the vicinity of two-dimensional networks of topolog- 
ical obstacles shows interesting behavior. We have investigated the diffusion constants of a 
polymer of polymerization index in a solution in which the two-dimensional networks are 
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stacked with uniform spacing d. 

Three regimes of dynamics have been shown to exist as a function of N and d. 

Reptation regime When d is much smaher than the radius of gyration of the polymer N^^'^b, 
the dominant mode of relaxation is reptation and simple application of de Gennes's rep- 
tation theory reproduces our Monte Carlo data quite well. The diffusion constant is 
proportional to N~'^ and to d. 

Intermediate regime While the chain is almost free in the direction parallel to the mem- 
brane when d > N^^'^b, the motion of the chain perpendicular to the membrane is still 
dominated by membrane penetration up to d ~ Nb. Hence the dynamics of the chain 
is highly anisotropic in this regime, the ratio D^y/Dz being R^Dji/ {dD-^) ~ N/d. 

Rouse regime When d becomes larger than the chain length Nh the diffusion of the chain 
is predominantly Rouse-like and the correction to the Rouse diffusion constant is 
0{M^^^TZ^/\) for and 0{TZ^/\) for I?xy The penetration time for a single isolated 
membrane is found to be proportional to N^^"^. This can be interpreted as slowing 
down of the diffusion of a polymer in the perpendicular direction to the membranes by 
an extra time of N'^^'^ti per membrane. 

This study in principle can provide a new way of filtrating polymers using highly perme- 
able membranes. The permeability for a single membrane has been found to be proportional 
to N^"^ . The permeability for multiple membranes will be roughly proportional to divided 
by the thickness of the stack, providing very rich behavior depending on the spacings between 
membranes. 

Even though we hope that this study gives impetus to filtration techniques using fishnet- 
like membranes, this work omits important effects. For example we have completely neglected 
excluded volume effects, which could, for instance, make the membranes effectively repulsive 



to the polymers. The effect of membrane fluctuations ||2^, 21, 22 1 were also neglected in this 
work. Even though the time scale of large length scale fluctuations of a two-dimensional 
system is much larger than that of linear polymer relaxation, it could change some important 
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characteristics of the filtration process, for instance by varying the distances between the 
membranes, which were assumed uniform in this paper. 
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Figure 1: The logarithm of penetration time, r, defined in the text is plotted against the 
logarithm of the polymerization index, N . The membranes are placed at distances 5, 6, 8, 
10, 12, 14 and 18 for AT = 8, 10, 16, 24, 36, 50 and 80. The distances are chosen so that the 
probability of the chain touching two membranes at the same time is extremely small. 
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Figure 2: The diffusion constants are given as functions of d for = 50. The circles represent 
Dxy and the triangles D^. The Rouse value is around 5.5 x 10~^ in our unit. Each pair of 
points represents an average over 1000 Monte Carlo runs. 
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Figure 3: Scaling graph for for N = 10(open triangle), 16(cross), 24(pentagon), 36(filled 
triangle), 50(snowflake) and 80(circle). The scaling function is = Di/N^-^F{d/Rg) where 
Rg = is used. F is pretty linear when d/Rg is small and the scaling breaks down 

beyond that. Note that each data point was obtained separately. Each point is an average 
of 1000 Monte Carlo runs. 
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Figure 4: Scaling graph for Dxy The symbols represent the same A^'s as in Fig. |. The 
break-down of the scaling when d > Rg is more dramatic here than in Fig. |3|. 
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Figure 5: Average mean square displacements of 40 monomers at the center of a chain with 
N = 200. Ti represents z-component (triangle) or x(or ?/)-component(circle) of the coordinate 
of a monomer. The curves were averaged over 200 Monte Carlo runs. After the t^l"^ regime 
-Dz seems to show t^/^ behavior again, but it is not very clear from this graph. 
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